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Abstract. We provide a new derivation of the conditions of dynamical and thermodynamical stability of 
homogeneous and inhomogeneous isothermal distributions in the Hamiltonian Mean Field (HMF) model. 
This proof completes the original thermodynamical approach of Inagaki [Prog. Theor. Phys. 90, 557 
(1993)]. Our formalism, based on variational principles, is simple and the method can be applied to more 
general situations. For example, it can be used to settle the dynamical stability of polytropic distributions 
with respect to the Vlasov equation [Chavanis & Campa. larXiv:1001.2109] . For isothermal distributions, 
the calculations can be performed fully analytically, providing therefore a clear illustration of the method. 

PACS. 5.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and chaos - 05.20.Dd Kinetic 
theory - 64.60.De Statistical mechanics of model systems 



1 Introduction 

The statistical mechanics of systems with long-range inter- 
actions has recently been the object of an intense activity 
|1I2I3I4| . A system with long-range interactions is charac- 
terized by a binary potential u{r) which decreases at large 
distances r slower than r~'^ where d is the dimension of 
space. Systems with long-range interactions are numerous 
in nature and include for example self-gravitating systems, 
two-dimensional vortices, bacterial populations experienc- 
ing chemotaxis, neutral and non neutral plasmas, wave- 
particle systems, free-electron lasers,... These systems ex- 
hibit a very interesting dynamics and thermodynamics. 
On a dynamical point of view, they display robust and 
long-lived quasi-stationary states (QSSs) that are non- 
Boltzmannian. These distributions are stable steady states 
of the Vlasov equation on the coarse-grained scale result- 
ing from a process of violent relaxation [S] . On a thermo- 
dynamical point of view, their Boltzmannian statistical 
equilibrium states can display numerous types of phase 
transitions due to ensembles inequivalence |6I7] . Ensem- 
bles inequivalence is generic for systems with long-range 
interactions, as first evidenced in astrophysics (see reviews 
in |8l9ll0j ). but not compulsory. 

A simple toy model of systems with long-range interac- 
tions, called the Hamiltonian Mean Field (HMF) model, 
has received a particular attention. This model consists 
of N particles of unit mass moving on a ring and inter- 



acting via a cosine potential u = N cos( 



j) where 



9i denotes the angle that makes particle i with an axis 
of reference. This model was first introduced by Messer 
& Spohn [TT] and called the cosine model. It was reintro- 



duced independently ten years later by different groups 
of researchers |12I13I14I15I16J and was very much studied 
since then (see a short history in [13 )• In particular, the 
seminal paper of Antoni & RufFo |16| inspired many works 
on the subject. 

The HMF model exhibits two successive types of re- 
laxation. For short timescales, the distribution function 
is governed by the Vlasov equation that ignores correla- 
tions (or "collisions") between particles. In this regime, 
the system can experience a collisionless violent relax- 
ation towards a steady state of the Vlasov equation on 
the coarse-grained scale. In principle, this quasistationary 
state (QSS) can be predicted by using Lynden-Bell's sta- 
tistical theory of violent relaxation 118119120121"!^ . How- 
ever, this prediction may fail because collisionless relax- 
ation can be incomplete [22123124] . Since the Vlasov equa- 
tion admits an infinite number of steady solutions, the 
prediction of the QSS actually reached by the system in 
case of incomplete relaxation is very difficult, or even im- 
possible. Nevertheless, it can be useful to derive general 
stability criteria in order to determine which distribu- 
tions are stable or unstable with respect to the Vlasov 
equation. This is the problem of Vlasov dynamical stabil- 
ity |l3ll.qilfill7l24l2ftl26l27l28[^ . Of course, only stable 
states are relevant to characterize QSSs in the context 
of violent relaxation. On a longer timescale, the system 
is expected to achieve a Boltzmannian statistical equi- 
librium state due to the development of correlations be- 
tween particles (finite N effects or graininess). This sta- 
tistical equilibrium state corresponds to the distribution 
that maximizes the Boltzmann entropy at fixed mass and 
energy. Only global or local entropy maxima are relevant 
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(mininia or saddle points must be discarded). This is the 
problem of thermodynamical stability |11I14I15I16I17I30] . 
These problems of dynamical and thermodynamical sta- 
bility have been investigated in the past using different 
methods that we shall briefly review in this introduction. 

Let us first discuss the problem of thermodynamical 
stability. Messer &; Spohn [TT] considered a potential en- 
ergy of the form U = N~'^ Si<j ^(^^i, Xj) and proved rig- 
orously that the mean field approximation is exact for 
N — > +CX) and that the statistical equilibrium state in the 
canonical ensemble corresponds to the global minimum 
of free energy F[f] at fixed mass M. Considering specif- 
ically the cosine model, they showed that it displays a 
second order phase transition from a homogeneous phase 
to a clustered phase below a critical temperature Tc. In- 
agaki [T3] studied the thermodynamical stability of the 
modified Konishi-Kaneko |12] system in the microcanoni- 
cal ensemble by considering the maximization of entropy 
S[f] at fixed mass M and energy E. By using the theory of 
Poincare on linear scries of equilibria (see, e.g., |10I31| for 
details) , or by studying the sign of the second order vari- 
ations of entropy (for the homogeneous phase) , he showed 
that the system displays a microcanonical second order 
phase transition at a critical energy Ec, corresponding to 
the critical temperature TJj- Similar results were obtained 
by Pichon [TS] who showed in addition that this type of 
phase transitions could explain the formation of bars in 
disk galaxies. Antoni & Ruffo [TB] studied the statistical 
equilibrium state of the HMF model in the canonical en- 
semble directly from the partition function. They simpli- 
fied it by using the Hubbard-Stratonovich transformation 
and the saddle point approximation valid for N — > -fcxD. 
They evidenced a second order phase transition at T = Tc 
in the canonical ensemble and performed numerical simu- 
lations in the microcanonical ensemble. These simulations 
show a discrepancy with the theoretical results close to 
the transition energy Ec, but this discrepancy is not due 
to ensembles inequivalence but to nonequilibrium effects 
|16I24I32] . More recently, Barre et al. |5D] studied the sta- 
tistical mechanics of the HMF model by applying large 
deviation technics. They confirmed the existence of mi- 
crocanonical and canonical second order phase transitions 
and the equivalence of statistical ensembles. Finally, Cha- 
vanis et al. |17| pursued the thermodynamical approach 
of Inagaki [14j based on variational principles. In partic- 
ular, they reduced the stability problem to the study of 
an eigenvalue equation and solved this eigenvalue equa- 
tion numerically for any energy and analytically close to 
the critical point {Ec,Tc). They proved by this method 
that the statistical ensembles are equivalent and that the 
homogeneous states are stable for E > Ec (or T > Tc) 
and unstable for E < Ec (or T < Tc). Below that critical 
energy or critical temperature, they are replaced by in- 
homogeneous states that arc always stable. This method 
has the advantage of showing which type of perturbation 



^ Although not explicitly mentioned by Inagaki [TJ, it is 
clear that a direct application of the Poincare theorem shows 
that the statistical ensembles are equivalent for this model. 



is able to trigger the instability of the homogeneous phase 
below Ec or Tc. 

Let us now review the results concerning the dynam- 
ical stability of steady states of the Vlasov equation in 
the context of the HMF model. We first consider the lin- 
ear dynamical stability problem. Inagaki & Konishi |13j 
and Pichon |15] studied the linear dynamical stability of 
the isothermal (Maxwell) distribution with respect to the 
Vlasov equation. They considered the homogeneous phase 
and derived the dispersion relation using the methods of 
plasma physics and stellar dynamics. They showed that 
the system becomes dynamically unstable below the crit- 
ical temperature Tc or the critical energy Ec (the same as 
the ones arising in the thermodynamical approach) lead- 
ing to an instability similar to the Jeans instability in self- 
gravitating systems. They showed that only the modes 
n ^ ±1 grow, leading to the formation of a single clus- 
ter. Inagaki & Konishi ^3 ^-l^o compared their theoretical 
results with direct numerical simulations, finding a good 
agreement. Antoni & Ruffo |16| studied the linear dynami- 
cal stability of the spatially homogeneous waterbag distri- 
bution and determined the critical temperature T^ and the 
critical energy E'c marking the separation between stable 
and unstable states. More recently, Choi & Choi [551 j Cha- 
vanis et al. [T7] and Jain et al. [57] completed these studies 
and obtained explicit expressions for the growth rate and 
pulsation period of isothermal, polytropic and waterbag 
distributions. On the other hand, Chavanis & Delfini [55] 
performed an exhaustive study of the linear dynamical 
stability of the HMF model by using the Nyquist method. 
They considered various types of distributions (single and 
double humped) and derived general stability criteria and 
stability diagrams for both attractive (ferromagnetic) and 
repulsive (antiferromagnetic) interactions. 

Let us finally review the results concerning the for- 
mal nonlinear dynamical stability of a steady state of the 
Vlasov equation. A distribution function is said to be for- 
mally stable |33j if it is a local minimum or a local maxi- 
mum of an energy-Casimir functional (i.e. the second vari- 
ations of the energy-Casimir functional are positive defi- 
nite or negative definite). Formal stability implies linear 
stability, but the converse is wrong in general. Yamaguchi 
et al. [25j derived a necessary and sufficient condition of 
formal nonlinear dynamical stability for spatially homo- 
geneous distribution functions of the HMF model. They 
observed that linear stability and formal stability criteria 
coincide in that case. Chavanis et al. |17| reconsidered the 
formal stability problem on a new angle (see also [28134] 
for more details) which can be extended to more general 
situation^ They showed that the variational problem for 
the distribution function f{9,v) is equivalent to a simpler 
variational problem for the density p{9). This equivalence 
is valid for both homogeneous and inhomogeneous dis- 
tributions. For spatially homogeneous distributions, they 
showed that the condition of formal stability can be writ- 
ten as a condition on the velocity of sound in the corre- 
sponding barotropic gas. This is equivalent to the criterion 



^ This method is related to the Antonov first law in astro- 
physics |35I36| . 
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derived by Yamaguchi et al. [2S] but expressed in a differ- 
ent manneiij. For spatially inliomogeneous distributions, 
tfiey reduced tlie formal stability problem to the study of 
an eigenvalue equation. Cliavanis & Delfini [55] showed 
that this eigenvalue equation can be solved analytically 
at the point of neutral stability, providing therefore an 
explicit condition to locate this point in the series of equi- 
libria. They also generalized the preceding results to arbi- 
trary potentials of interaction u(r,r') and discussed more 
refined criteria of dynamical stability that are obtained by 
conserving a larger class of Vlasov invariants. 

Recently, Campa & Chavanis [52] derived a very gen- 
eral criterion of linear dynamical stability valid for homo- 
geneous and inhomogeneous distributions. They also de- 
rived sufficient conditions of linear and formal dynamical 
stability that arc weaker than the general criterion but 
more explicit. For spatially homogeneous distributions, 
they proved that the criteria of formal and linear stabil- 
ity coincide. As a by-product, their results also return the 
thermodynamical stability criteria obtained previously by 
different methods. 

Let us finally mention that the Vlasov dynamical sta- 
bility of polytropic distributions has been studied by Cha- 
vanis & Campa [24j by plotting the series of equilibria 
and using the Poincare criterion. On the other hand, the 
Vlasov dynamical stability and the thermodynamical sta- 
bility of the Lynden-Bell distributions have been studied 
by Antoniazzi et al. |20j and Staniscia et al. (22j by solving 
the variational problem numerically. These studies exhibit 
a rich phase diagram with several types of phase transi- 
tions showing the complexity of the stability problem in 
the general case. 

In this paper, we shall present a new method to de- 
termine the dynamical and thermodynamical stability of 
isothermal distributions of the HMF model. The idea is 
to start from the fundamental variational problems (jl2[) 
and (J13p and transform them into equivalent but simpler 
variational problems until a point at which they can be 
explicitly solved. This completes the thermodynamical ap- 
proach of Inagaki [Tl] by proving analytically the stability 
of the inhomogeneous phase which was not done in Ina- 
gaki's paper. An interest of our approach is its simplicity 
and generality. Indeed, it can be extended to solve Vlasov 
and Lynden-Bell stability problems. For example, it has 
been used recently to settle the Vlasov dynamical stabil- 
ity of polytropic distributions (see Sec. 8 of [H]). How- 
ever, in that case, the calculations are less explicit than 
for isothermal distributions. It is therefore interesting to 
develop the calculations in detail in the case of isothermal 
distributions (where they are fully analytical) in order to 
clearly illustrate the method. Although we rederive known 
results in a different manner, we think that the present ap- 
proach is interesting and potentially useful to tackle more 
general problems. In addition, our approach not only de- 
termines the strict equilibrium state (global maximum of 



^ These criteria have been used in |17l25j to determine the 
formal stability of isothermal and polytropic distributions and 
in [18] to determine the formal stability of the Lynden-Bell (or 
Fermi-Dirac) distribution. 



entropy) but it can also say whether a critical point of en- 
tropy is metastable (local entropy maximum) or unstable 
(saddle point of entropy). Although isothermal distribu- 
tions of the HMF model do not display metastable states, 
this information can be valuable in more general situations 
|20l22l24j where our method can be applied. 



2 Series of equilibria 

The HMF model is a system of N particles of unit mass 
m = 1 moving on a circle and interacting via a cosine 
potential. The dynamics of these particles is governed by 
the Hamilton equations |14I15I16I17J : 



dt 

N 



dH 

dvi' 



dt 



H 



1 ^ — ^ 2 '^ \ — ^ 



dH 



(1) 



j#j 



where 6i e [— tt, tt] and — oo < Vi < -l-cio denote the 
position (angle) and the velocity of particle i and k is 
the coupUng constant (we assume here that fc > 0). 
The thermodynamic limit corresponds to W — > +oo in 
such a way that the rescaled energy e — SnE/kM^ re- 
mains of order unity. We can take k ~ 1/A^ which is the 
Kac prescription. In that case, the energy is extensive, 
E/N ^ 1, but non-additive. For N -^ +oo, the mean 
field approximation is exact and the A'^-body distribu- 
tion function is a product of N one body distributions: 

PN{0i,Vi,...,9N,VN,t) = Pi{0i,Vi,t)...Pi{9N,VN,t). 

Let us introduce the distribution function f(9,v,t) = 
NPi{d, V, t). For a fixed interval of time and N — >■ -|-oo, the 
evolution of the distribution function f{9, v, t) is governed 
by the Vlasov equation 



dl dl_d^dl^ 
dt 89 39 dv 



where 



k 

2^ 



27r 



^i^,t) = -7r cos{e -6' )p{9',t)de', 



(2) 



(3) 



is the self-consistent potential generated by the density of 
particles p{9, t) = J f{9, v, t) dv. The mean force acting on 
a particle located in 9 is (F) = —d<P/d9{9,t). Expanding 
the cosine function in equation ([3]), we obtain 



wher 



<?(6», t) = B^ cos 9 + By sin 6*, 



fc f 
B^. = -—- / pie, t) cos 9 d9, 
2tt J 



By = -— / p{9,t) sine d9, 



(4) 



(5) 



(6) 
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are proportional to the magnetization (with the opposite 
sign). The magnetization can be viewed as the order pa- 
rameter of the HMF model. 
Let us introduce the mass 



M[p]^ J pd9, 



(7) 



and the mean field energy 



/ fiP dodv + - I p<Pde = K + W, (8) 



E[f]^^ I fv^dOdv . 2 



where K is the kinetic energy and W the potential en- 
ergy. Using equations (|4])-(|6l), the potential energy can be 
expressed in terms of the magnetization as 



W 



7rB2 
k 



We also introduce the Boltzmann entropy 



S[f] 



/In(^) d0dv, 



and the Boltzmann free energy 



F[f]^E[f]-TS[f], 



(9) 



(10) 



(11) 



where T = 1/(3 is the temperature. In the microcanonical 
ensemble, the statistical equilibrium state is determined 
by the maximization problem 



max{S[f]\E[f] 



E, M[f] - M} . 



(12) 



In the canonical ensemble, the statistical equilibrium state 
is determined by the minimization problem 



mm{F[/]|M[/] = M}, 



(13) 



The Boltzmann entropy functional pop and the maximum 
entropy principle (|12p can be justified from a standard 
combinatorial analysis (see, e.g., |10ll7j V The distribution 
function f{9, v) that is solution of p2p is the most probable 
macroscopic state, i.e. the macrostate that is the most 
represented at the microscopic level, assuming that the 
accessible microstatcs (with the proper values of mass and 
energy) are equiprobable (see Appendix ICj) . 

We shall first determine the critical points of these vari- 
ational problems. This will allow us to set the notations 
that will be needed in the following. The critical points 
of the maximization problem (J12p are determined by the 
variational principle 



SS - (3SE - aSM = 0, 



(14) 



where /3 = 1/T and a are Lagrange multipliers associated 
with the conservation of energy and mass. The critical 
points of the minimization problem (J13p are determined 
by the variational principle 



5F + aT5M = 0, 



(15) 



where a is a Lagrange multiplier associated with the con- 
servation of mass. Since T is fixed in the canonical ensem- 
ble, it is clear that equation (jlSp is equivalent to equation 
(fn)) . Therefore, the optimization problems ([T^ and ([T5)) 
have the sam,e critical points. Performing the variations in 
equations ([T^ and p^ , we find that the critical points are 
given by the mean field Maxwell-Boltzmann distribution 



fiO,v)^A'e 



-f)[4+${e)] 



(16) 



where A' = e ^ " is a constant. Integrating over the ve- 
locity, we get the mean field Boltzmann distribution 



p(0)=Ae-'^*W, 



(17) 



where A ~ {2tt/I3)^/^A' . Using the expression ^ of the 
potential, the distribution function p6l) can be rewritten 



fie,v)^A'e 



-Pl^+Bj: cose+S„ sine) 



It is convenient to write B^ 



Bcosd) and B,, 



(18) 



i?sin( 



with B = {Bl 
pression takes the form 



^2)1/2 



In that case, the foregoing ex- 



f{0,v) 



yi'g-'3[4+-Bcos(e-0)] 



The corresponding density profile is 



(19) 



(20) 



The amplitude A and the magnetization B are determined 
by substituting equation (pn|) in equations ([5]), ([5]) and ([7]). 
This yields 

A^^-^L^, (21) 



and 



2^/o(/3B)' 

2-kB _ Ii{l3B) 
Jm ~ loWB)' 



(22) 



where /„(x) is the modified Bessel function of order n. 
Equation (|22p determines the magnetization i? as a func- 
tion of the temperature T. Then, A is given by equation 
((2T|) . Note that the critical points are degenerate. There 
exists an infinity of critical points which differ only by 
their phase </>, i.e. by the position of the maximum of the 
density profile. They have the same value of entropy or 
free energy (see below). In the following, we shall take 
(p = without loss of generality. In that case, B^ ^ B and 
By ~ 0. Then, the distribution function and the density 
can be written 



m-) = Uz 



1/2 



Pi0)e' 



p{0) 



M 



-pBcose 



2TrIoipB) 



(23) 



(24) 



where B is determined in terms of T by equation 
The study of the self-consistency relation ([2^ is classical: 



P.H. Chavanis: Dynamical and thermodynamical stability of isothermal distributions in the HMF model 



i? = is always solution while solutions with B ^ {) only 
exist for T <Tc = kM/Air J4I14I17J . 

Let us now determine the expressions of the energy, 
entropy and free energy. For the Maxwell-Boltzmann dis- 
tribution (|23p . the kinetic energy is 



3f 



1 



-MT. 



(25) 



Combining this relation with equation 0, we find that 
the total energy E = K + W is given by 



E ^ -MT 
2 



ttB^ 







1 


- 


-1 


B=0 




- 


- I 


(unstable) 




. 


\ B>0 








- \ (stable) 


1 






- v_ 


\ E,>n, 


B=0 


- 


" e , 




(stable) 




1 





(26) 



1 2 

e=87iE/kM' 



The series of equilibria giving T as a function of E is Fig- 1- Series of equilibria (caloric curve) giving the tempera- 
determined by equations i^ and dH]) by eliminating B. ture as a function of energy. 
The relation that gives the magnetization _B as a function 
of the energy E is determined by equations ([^ and (PSJ) 

by eliminating T. Finally, using equations (UHl) and dMD, These rclatfons apply to the inhomogcneous states (b^O). 

the entropy is given by Po,. ^^e homogeneous states (6 = 0), we have 



1 



5*= -MlnT 
2 



plnpdO, 



(27) 



1 



■ In 77. 



(35) 



77 2 V V 

up toa term ^M + iMln(2^) + MlnM. Using equation r^^^ magnetization b takes values between and 1, the 
((231), it can be rewritten inverse temperature 77 between and -hoo and the energy 

., r, n2 e between emin — —2 and -l-oo. The homogeneous states 

" '" -'""" ^^23^ exist for any 77 > and for any e > 0. The inhomogcneous 

states exist for any t] > r]c = I and any emin ^ e ^ Cc = 1- 



S=-MlnT + MlnIoiPB)- ^^ 



up to a term ^M + |A/ ln(27r). The relatfon between the The bifurcatfon point is focated at 

entropy S and the energy E is determined by equations 

(Pg)) . (pS)) and dm. Using equations dH]) and dig), the , 

free energy (|lip is given by 



SnEc 
kA'P 



= 1, 



F = -MT - -MTlnT- MT In Iq{I3B) 



ttB' 



Close to the bifurcation point (7/ 
(29) 1~), we get (see Appendix IXj) : 



l3ckM 
An 

-^ Vc = 



= 1. 

= 1+, 



(36) 



up to a term -iMT- |Mnn(27r). The relation between ^^ /2(?7 - 1) 6' 
the free energy F and the temperature T is determined 
by equations dH]) and dH])- 

It is convenient to write these equations in parametric 
form by introducing the parameter x = /3B. Then, we 
have 



5(1 -e), V 



h^ 



b = 



2ttB _ h{x) 
Jm ~ /o(.t)' 

PkM _ X 
47r ~ 2b{x) ' 



sc^--il-e), /^l--(77-l)^ 

^ ' Close to the ground state (77 — >> -t-00, e -^ emin 
get (see Appendix [X|) : 



(31) ^^1-^, b' 



(37) 

(38) 
-2), we 

(39) 



87r£; 
k^P 



1 



77(2;) 



26(x•)^ 



(32) 



s~ln(e + 2), /: 



■ In 77. 



(40) 



^ = --hi77(a;) +ln/o(a;) - 2j^{x)b{xf , 



8nF 
kM^ 



e{x) 



ri{x) 



s{x). 



(33) 



(34) 



From these relations, we can obtain the curves T{E), 
B{T), B{E), S{E) and F{T) in parametric form. These 
curves are plotted in Figures [1][S] for completeness. These 
results are well-known and they have been derived in 
many papers [4114116117130] using different methods. The 
present approach, that complements the original approach 
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2 3 

T|=|3kM/47t 

Fig. 2. Magnetization (order parameter) as a function of tem- 
perature. 




2 3 

ri=|3kM/43t 

Fig. 5. Free energy as a function of temperature. 




Fig. 3. Mai 

ergy. 



e=8JiE/kM 
gnetization (order parameter) as a function of en- 




e=83tE/kM 
Fig. 4. Entropy as a function of energy. 



of Inagaki [T3] , is the most direct and the most complete. 
Indeed, these relations characterize all the critical points 
of p^ and ([T^ . Now, we must select among these critical 
points those that are (local) maxima of S at fixed E and 
M (microcanonical ensemble) and those that are (local) 
minima of F at fixed M (canonical ensemble). This is the 
object of the next sections. 



Remark: using the Poincare theorem (see, e.g. |10I31| ). 
we can directly conclude from the series of equilibria 
that the homogeneous states are microcanonically (resp. 
canonically) stable for E > Ec (resp. T > Tc) while 
they are microcanonically (resp. canonically) unstable for 
E < Ec (resp. T < Tc). On the other hand, the inhomoge- 
neous states are always stable. The caloric curve therefore 
exhibits a microcanonical second order phase transition 
marked by the discontinuity of f3'{E) = S"{E) &i E = Ec 
and a canonical second order phase transition marked by 
the discontinuity of E'{f3) = (l3F)"{fi) at /3 = ^c- The en- 
sembles are equivalent. In order to illustrate our method, 
which remains valid in more general situations, we shall 
however treat both canonical and microcanonical ensem- 
bles. 



3 Canonical ensemble 

3.1 The functionals F[f] and F[p\ 



The minimization problem (|13p has several interpreta- 
tions: 

(i) It determines the statistical equilibrium state of 
the HMF model in the canonical ensemble. In that ther- 
modynaniical interpretation 5* is the Boltzmann entropy, 
F is the Boltzmann free energy and T is the thermody- 
naniical temperature. The minimization problem (|13p can 
therefore be interpreted as a criterion of canonical ther- 
modynaniical stability. For isolated systems that evolve at 
fixed energy, like the HMF model, the canonical ensemble 
is not physically justified. For such systems, the proper 
ensemble to consider is the microcanonical ensemble and 
the statistical equilibrium state is given by (|12p . However, 
the canonical ensemble always provides a sujficient con- 
dition of microcanonical stability [6] and it can be useful 
in that respeclQ. Indeed, if a system is canonically stable 



* For systems with short-range interactions, the ensembles 
are equivalent and we can choose the one that is the most 
convenient to make the calculations. For systems with long- 
range interactions, the ensembles may not be equivalent: grand 
canonical stability implies canonical stability which itself im- 
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at the temperature T, then it is automatically granted to 
be niicrocanonically stable at the corresponding energy 
E ~ E{T). Therefore, we can start by this ensemble and 
consider the microcanonical ensemble only if the canoni- 
cal ensemble does not cover the whole range of energies. 
On the other hand, for systems in contact with a thermal 
bath fixing the temperature, like the Brownian Mean Field 
(BMF) model [T7|, the canonical ensemble is the proper 
ensemble to consider and the statistical equilibrium state 
is given by (fT3|) . 

(ii) It determines a particular steady state of the 
Vlasov equation that is formally nonlinearly dynamically 
stablqj- The minimization problem (jl3|) can therefore be 
interpreted as a sufficient condition of dynamical stabil- 
ity. In that dynamical interpretation, 5 is a particular 
Casimir (pseudo entropy), F is an energy-Casimir func- 
tional (pseudo free energy) and T is a positive constant. 
It can be convenient to develop a thermodynamical anal- 
ogy |34] to study this dynamical stability problem and use 
a common vocabulary. In this way, the methods devel- 
oped in thermodynamics can be applied to the dynamical 
stability context. 

It is shown in Appendix A. 2. of [37] that the solution 
of (fT3|) is given by 



f{0,v)^ 



1/2 



p{B)e- 



where p{0) is the solution of 

mhY{F[p]\M[p\ = M}, 



F[p\ = \ I p'^dB + T I plnpdB. 



(41) 



(42) 



(43) 



where 



Therefore, the minimization problems (|13|) and (|42|) are 
equivalent: 

131)^ dH. 



(44) 

This equivalence holds for global and local minimization 
[37] : (i) f{B, v) is the global minimum of ^ iff p{e) is the 
global minimum of (|^^ and (ii) f{6, v) is a local minimum 
of (fT5|) iff p{B) is a local minimum of (02]). We are therefore 
led to considering the minimization problem (j42[) which is 
simpler to study since it involves the density p{B) instead 
of the distribution function f(v,B). 

Before that, let us compare the conditions of stabil- 
ity issued from ((T^ and P^ . The critical points of (|13l) 
are given by equations ^^ , ([M)) and (f^ where /? is pre- 
scribed. A critical point of F[f] at fixed mass is a (local) 
minimum iff 



^F= i f 6pS^dB + T I ^-^dBdv> 



0, 



(45) 



plies microcanonical stability, but the converse is wrong in gen- 
eral. 

^ This is a refined condition of formal stability with respect 
to the usual criterion [33] since the mass is treated here as a 
constraint (see [28I29J for a more detailed discussion). 



for all perturbations 5/ that do not change the mass: 
5M = 0. On the other hand, the critical points of (|42]) are 
given by equations 1^^ and (|22p where j3 is prescribed. A 
critical point of F[p] at fixed mass is a (local) minimum 
iff 



S'F^ 



1 



6pS<P dB + T 



M 

2p 



d6>Q, 



(46) 



for all perturbations 5p that conserve mass: 5M — 0. This 
stability criterion is equivalent to the stability criterion 
(|45p but it is simpler because it is expressed in terms of 
the density instead of the distribution function |37| . 

Remark 1: the thermodynamical approach of Messer & 
Spohn [TT] in the canonical ensemble directly leads to the 
minimization problem (|42p for the density, and justifies it 
rigorously. Our approach recovers this result by another 
method. It also shows that this minimization problem pro- 
vides a sufficient condition of dynamical stability with re- 
spect to the Vlasov equation. 

Remark 2: considering the BMF model [17) . the mini- 
mization problem (|T3| determines stable steady states of 
the mean field Kramers equation and the minimization 
problem (|42p determines stable steady states of the mean 
field Smoluchowski equation. According to the equivalence 
(jii]) . we conclude that a distribution function f{9,v) is a 
stable steady state of the mean field Kramers equation 
iff the corresponding density field p{9) is a stable steady 
state of the mean field Smoluchowski equation (see [55] 
for a more general statement). 

Remark 3: the equivalence between ([13]) and ([42]) can 
be extended to a larger class of functionals of the form 
S[f] = — J C{f) dBdv where C is convex. Such functionals 
(Casimirs) arise in the Vlasov dynamical stability problem 
[28I29I33J . We refer to [17134139] for a detailed discussion 
of this equivalence. 



3.2 The function F{B) 

3.2.1 Global minimization 

The equivalence ([^^ is valid for an arbitrary potential of 
interaction u{r, r'). Now, for the HMF model, the problem 
can be further simplified. Indeed, the potential energy is 
given by equation ^ so that the free energy ([^ can be 
rewritten 



F[p] 



ttB^ 
k 



T / phypdB. 



(47) 



We shall first determine the global minimum of free en- 
ergy at fixed mass. To that purpose, we shall reduce the 
minimization problem (|47p to an equivalent but simpler 
minimization problem. 

To solve the minimization problem (|47p . we proceed in 
two stepCl- we first minimize F[p\ at fixed M and B^ and 



® We have used this method in different situations (see, e.g. 
134136137) 1. 
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By. Writing the variational principle as 



plnpde] + aSM + pJB^ + p.y5By = 0, (48) 



we obtain 



,-A^ 



3e-A„ 



where A 



e 



Pi{0)=Ae 

A 
s I 
jy. If we write Bx ~ B cos 
fir 
A = (A^ + Ay)^/^. Equation (|49]) can be rewritten 



(49) 
The 



IfM:^ and Xy = -^fiy 

Lagrange multipliers are determined by the constraints M, 
Bx and By . If we write B^ ~ B cos cf) and By ~ B sin 
then we find that Aa; = A cos and Ay = A sin (f> with 



(50) 



Finally, A and A are determined in terms of M and B 
through the equations 



A = 



M 



277 Jo (A) ^ 



and 



b = 



2ttB _ /i(A) 



(51) 



(52) 



kM Jo (A)' 

Equation ([501 is the (unique) global minimum of F[/9] with 
the previous constraints since S'^F = ^T J ^ ^' d9 > 
(the constraints are linear in p so that their second varia- 
tions vanish). Then, we can express the free energy F[p\ as 
a function of B by writing F{B) = F[pi]. After straight- 
forward calculations, we obtain 



F{B) 



nB' 
k 



T\ 



2ttB 



MT\nIo{X), (53) 



where X{B) is given by equation (|52p . Finally, the mini- 
mization problem (j42p is equivalent to the minimization 
problem 



min {F{B)}., 



(54) 



in the sense that the solution of P^ is given by equations 
dSni), dSH and dSH) where B is the solution of dM]). Note 
that the mass constraint is taken into account implicitly 
in the variational problem ([M]) . Therefore, (|^^ and (|54p 
are equivalent for global minimization. However, (|54p is 
much simpler because, for given T and M , we just have 
to determine the minimum of a function F{B) instead of 
the minimum of a functional F[p] at fixed mass. 

Let us therefore study the function F{B) defined by 
equations ([55]) and (|52p . Its first derivative is 



F'{B) 



2TrB 



TX 



2n 



MT 



(^ttB /^(A) 
\kM Jo (A) 



dX 
dB' 



(55) 



Using the identity /o(A) = /i(A) and the relation ((52)) . we 
see that the term in parenthesis vanishes. Then, we get 



F'{B) 



27r 



{TX - B) 



(56) 




Fig. 6. Graphical construction determining the critical points 
of F[B) and their stability. The critical points are determined 
by the intersection(s) between the curve h — 6(A) defined by 
equation (|52|) and the straight line b — \/{2r}). The critical 
point is a minimum (resp. maximum) of F{B) if the slope of 
the curve 6(A) at that point is smaller (resp. larger) than the 
straight line b = \/{2rj). 



The critical points of F{B), satisfying F'{B) = 0, corre- 
spond therefore to 



X = x = PB. 

Substituting this result in equation ([5 
self-consistency relation 

2tiB _ hiPB) 



(57) 



we obtain the 



(58) 



which determines the magnetization _B as a function of 
the temperature T. This returns the results of Sec. [2] 

Now, a critical point oi F{B) is a minimum if F"{B) > 
and a maximum if F"{B) < 0. Differentiating equation 
([55)) with respect to B, we find that 



27r 



fix 



P'\B) = ^ ( ^^ - 1 



Therefore, a critical point is a minimum if 

dB 

<T 

dX ' 



(59) 



(60) 



and a maximum when the inequality is reversed. 

We can determine the minimum of the function F(B) 
by a simple graphical construction. To that purpose, 
we plot b = 2ttB /kM as a function of A according to 
equation ([5^ . This is represented in Figure [5] We find 
that 2tt B/kM — > 1 for A — >■ -l-oo. On the other hand, 
2ttB /kM ~ A/2 for A ^ 0. Therefore, the magnetization 
B takes values between and B^ax = -^^ According 
to equation ([57]) . the critical points of F{B) are deter- 
mined by the intersection of this curve with the straight 
line 2'KB/kM = {2TTT/kM)X. For T > T^ = ^, there is a 
unique solution B ^ corresponding to the homogeneous 
phase. For T < Tc, there are two solutions: a homogeneous 
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0.4 0.6 

2jiB/kM 

Fig. 7. Free energy F{B) as a function of the magnetization 
for a given value of the temperature. For T > Tc, this curve 
has a (unique) global minimum at B = 0. For T < Tc, this 
curve has a local maximum at B = and a global minimum 
at B{T) > 0. 



solution B ^ and an inhoniogcncous solution B(T) ^ 0. 
According to inequality (|60p . a solution is a minimum of 
F{B) if d{2TTB/kM)/dX < 2-KT/kM and a maximum if 
d{2wB/kM)/dX > 2'KT/kM. Therefore, a critical point of 
free energy F(B) is a minimum (resp. maximum) if the 
slope of the main curve is lower (resp. higher) than the 
slope of the straight line at the point of intersection. From 
this criterion, we easily conclude that: for T > Tc, the ho- 
mogeneous solution B = is the global minimum oi F{B); 
for T < Tc, the inhomogeneous solution B{T) 7^ is the 
global minimum of F{B) while the homogeneous solution 
B = is a local maximum. Now, using the equivalence 
between (|42)) and ((54)) for global minimization, we con- 
clude that the global minimum of the functional F[p\ at 
fixed mass is the homogeneous solution for T > Tc and the 
inhomogeneous solution for T < Tc. Furthermore, we will 
show in the next section that the homogeneous solutions 
for T < Tc are saddle points of F[p] at fixed mass. 

To complete our analysis, it can be useful to plot 
the function F{B) for a prescribed temperature. Using 
equations ((53|) and (|52|) . the normalized free energy / = 
ATiF/{k]\'P) can be expressed in terms of A according to 



/(A) = - 






Af§-ln/o(A) 



(61) 



Eliminating A between the expressions ([CTjl and (j52[) . we 
obtain the free energy f{b) as a function of the magne- 
tization b for a fixed value of the temperature T (more 
precisely, for given 77, these equations determine f{b) in a 
parametric form). For T > Tc and T < Tc, this function 
displays the two behaviors described above, as illustrated 
in Figure [3 For T ^ Tc so that A, 5 ^ 0, we find that the 
free energy takes the approximate form 



fib) 



--l)b' + ^b\ (,y 

rj J 4:7] 



(62) 



For rj > j]c = 1, we explicitly check that the minimum 
satisfying /'(6) = and f"{b) > is given by the first 
relation in equation p7p . 

Remark 1: In Appendix lEl we plot the second varia- 
tions of free energy F"{B{T)) (related to the variance of 
the magnetization) as a function of the temperature and 
recover the previous conditions of stability. 

Remark 2: since the phase <j) does not appear in the 
function (|53p . this means that the inhomogeneous minima 
of F{B) for T < Tc arc degenerate: there exists an infinity 
of minima that only differ in their phases (j). 

Remark 3: since canonical stability implies micro- 
canonical stability [6], and since the series of equilibria 
P{E) is monotonic (see Figure [T]), we conclude that the 
maximum of entropy at fixed mass and energy is the ho- 
mogeneous solution for E > Ec and the inhomogeneous 
solution for E < Ec. Since we cover all the accessible range 
of energies, we conclude that the ensembles are equiva- 
lent. We shall, however, treat the microcanonical ensem- 
ble specifically in Sec. U since the method can be useful 
in other contexts where the ensembles are not equivalent 
(see, e.g. [M]). 



3.2.2 Local minimization 

We shall now show that the minimization problems (|42p 
and (j54p are also equivalent for local minimization. To 
that purpose, we shall relate the second order variations 
of F[p] to the second variations oi F{B) by using a suitable 
decompositiorQ. 

A critical point of (|42p is determined by the variational 
principle 



SF + aTSM = 0, 



(63) 



where a is a Lagrange multiplier accounting for the con- 
servation of mass. This leads to the distribution 



M 



P = 



2TrIo{l3B) 



-PBcost 



(64) 



where we have taken (f) = without loss of generality 
(see Sec. H]). The magnetization B is obtained by substi- 
tuting equation ([M|) in equation (O leading to the self- 
consistency relation ((22|) . Using equations (|46)) and ([3]), 
this critical point is a (local) minimum of F at fixed mass 
iff 

S^F ^ ^^iiSB^r + my?) + \t j^^P^d6>0, 

(65) 

for all perturbations 5p that conserve mass: J 6pd9 = 0. 
The corresponding variations of magnetization are 



k 
5Bx = —7^ I Sp cos 9 d9, 



(66) 



^ This method was previously used in the context of the sta- 
tistical mechanics of the 2D Euler [40141) and Vlasov [34] equa- 
tions. 
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5By = -— / 5psuiede. 



We can always write the perturbations in the form 

Sp = Sp\i+Sp±, (68) 



(•gy-j For A =: /3B, according to the self-consistency relation 
((22)) . we have B{l3B) = B. Differentiating equation ^ 
with respect to A, we get 



2^ ovx^ /((A)/o(A)-/^(A)/i(A) 



where (5p|| =(/.* + 1^:1; cos + j/j^ sin 6*)^ and 5p^ = (5p-(5/9||. . , ., ,,.^ _ ,.. ,,.. _ i/r q\ 

The second condition ensures that all the perturbations ° o\ ) ^ i\ )' i\ ) ~ 2^ o\ ) 

are considered. We shall now choose the constants ^, Vx 



and Vy such that 



/2(A)) and /o(x) — hix) = j/i(x), and recalling equation 
((79)) . the foregoing relation can be rewritten 



(5p|| d6l = 0, 



(5i?^ = / (5p|| cos 6 d6, 



5B.U 



This implies that 



27r 



5p|| sin 6'(i6'. 



5px_ dO = 0, 



(69) 



(70) 



(71) 



(72) 



kM ^ ' A/o(A) k-^NP ' 



(81) 



Taking A = I3B and introducing the function X{B), which 
is the inverse of B{X), we obtain the identity 



27r 1 



l-^^^-^S. (82) 



kM\'(B) ^ /3B Io{l3B) BAP' 
Comparing this relation with equation (j78p , we obtain 
27r Att'^B'^ 



I = 



k\'{B) 



k^M 



(83) 



Solving equations ([71)1 . (|75|) and (|75|) for p, u^ and i/j,, and 
using the result ([55]) . wc find that 



/• 



(5pi cos edd = 0, I 5p±_ sin OdO^Q. (73) 

The conditions (p5)) . ([70)1 and ([7T|) lead to the relations 

(74) 



z^, = -X{B)5B,, 



2nBX'{B) ^^ 

M = nr: dBx 

^ kM 



2tt 
pM^—y.,B^Q, 



kVrr 

6B, =pB- -^I, 
Ztt 



SBy^^'^iM-Il 

ZTT 



where we have defined 



pcos 9d9. 



(75) 
(76) 

(77) 



27r SBy 

''^ ~ '~kJT~l- 

Therefore, the perturbation 6p\\ takes the form 
'27rB 



(84) 
(85) 



<5p|| = -X'{B) 



kM 
2n 



cos 61 ) p{e)6Bx 
SB,, 



k ^^ ' M-I 



(86) 



By construction. (5p|| and 5p±^ are orthogonal for the scalar 
product weighted by 1/p in the sense that 



6p\\5pi 



d9 = 0. 



(87) 



This forms a system of three algebraic equations that de- 
termines the three constants p, Vx and Vy. Using the equi- Indeed, we have 
librium distribution (|64p , we find after simple algebra that 

/ -^^IL^ de= {p + Vx cos 6 + Vy sin 9)5 p±_ d9 = 0, 
(78) -^ ^ -^ 



I = M 



1 him 



PB loWB) 



where we have used the identity /o(a;) — hix) — -Ii{x). where we have used equations (|7^ and (|75)) to get the last 
Let us define the function B{X) by the relation equality. As a result, we obtain 



27rB(A) /i(A) 



kM 



4(A)- 



(79) 



(Sp)' 



d9 



i^PW 



■d9 



{SpA 



d9. 



(89) 
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Using equations ([55| and ([55)1 . wc obtain after simplifica- 
tion 

(90) 

Therefore, the second order variations of free energy given 
by equation (|65p can be written 



S^F^-{TX'{B)-l)iSB,f 
k 



■K /27r T 



k\k M 
Using identity ([5^ . we obtain 



^-l)(5S.)^ + ^T|i^d^. (91) 



b''F^-F"{B){6B,f 



IT /27r T 



k\k M 



'--l)iSB,f^\TJ 



1^ /■ ihi:? 



de. (92) 



For the inhomogeneous phase i? 7^ 0, using equation (|78p 
and the self-consistency relation ([2^ . we get 



fc 



(93) 



In that case, equation (|M|) reduces to 

S^~F=-F"{B)iSB,)^ + ^T f ^^^de. (94) 

On the other hand, for the homogeneous phase 5 = 0, 
equation ([75)) leads to 



2 



Using equation ([206]), equation ([92]) reduces to 
J^^ = iF"(0) [(55.)^ + (55,)^] + ir| (^ 



d6i. 



(96) 



These relations show that p is a local minimum of F[p] at 
fixed mass iff i? is a local minimum of F{B). Indeed, if 
F"{B) > then S^F > since the last term is positive. 
On the other hand, if F"{B) < it suffices to consider 
a perturbation of the form ([68p with Sp± = and 6p\\ 
given by equation ([55]) to conclude that S^F < for this 
perturbation. This implies that p is not a local minimum 
of F[p\ since there exists a particular perturbation that 
decreases the free energy. This is the case for the homoge- 
neous solutions when T < Tc since they are local maxima 
of F{B). Therefore, ([^^ and ([M]) are equivalent for local 
minimization. Combining all our results, we conclude that 
the variational problems ([T^ . ([^^ and ([M]) are equivalent 
for local and global minimization: 



m^mii^iM 



(97) 



4 Microcanonical ensemble 

4.1 The functionals S[f] and ^ip] 

The maximization problem (1121) has several interpreta- 
tions: 

(i) It determines the statistical equilibrium state of the 
HMF model in the microcanonical ensemble. In that ther- 
modynamical interpretation S is the Boltzmann entropy. 
The maximization problem (|12p can therefore be inter- 
preted as a criterion of microcanonical thermodynamical 
stability. 

(ii) It determines a particular steady state of the 
Vlasov equation that is formally nonlinearly dynamically 
stablqj. The maximization problem ([T^ can therefore be 
interpreted as a sufficient condition of dynamical stability. 
In that dynamical interpretation, S" is a particular Casimir 
(pseudo entropy). As explained previously, it is convenient 
to develop a thermodynamical analogy [34j to study this 
dynamical stability problem and use a common vocabu- 
lary. 

It is shown in Appendix A.l. of [37] that the solution 
of ([T^ is given by 



f{0,v) 



1/2 



pie)e- 



(98) 



where /3 = l/T is determined by the energy constraint 



£; = -MT + W, 



and p{6) is the solution of 



(95) where 



m&x{S[p]\M[p] = A'/}, 
p 



S[p] = -M\nT - I plnpdd. 



(99) 



(100) 



(101) 



Eliminating the temperature thanks to the constraint 
([99p , wc can write the entropy in terms of p alone as 

S[p] = - I p\npde + ^M\n{E^W[p]). (102) 

Therefore, the maximization problems ([T2]) and (llOOp are 
equivalent: 



Wi^mm. 



(103) 



This equivalence holds for global and local maximization 
[57] : (i) f{9,v) is the global maximum of ^ iff p{9) is 
the global maximum of pOOp and (ii) f{9, v) is a local 
maximum of ([12]) iff p{d) is a local maximum of (|100p . We 
are led therefore to considering the maximization problem 



This is a refined condition of formal stability with respect 
to the usual criterion [H^ since both the mass and the energy 
are treated here as constraints (see [28I29| for a more detailed 
discussion) . 
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(jlOOp which is simpler to study since it involves the density 
p{9) instead of the distribution function f{v,9). 

Before that, let us compare the conditions of stability 
issued from ((12]) and (jlOOp . The critical points of (HI]) 
are given by equations dMl, dMl, (El]) and (gB]) where E 
is prescribed. Furthermore, a critical point of S at fixed 
mass and energy is a (local) maximum iff (see Appendix 
E]): 



S^S 



mi 

2/ 



d0dv -^f3 6p5<Pde < 0, (104) 



for all perturbations 6f that do not change the mass and 
the energy at first order: SM = SE = 0. On the other 
hand, the critical points of (jlOOp are given by equations 
(|M|) . (im and (P5)) where E is prescribed (see Appendix 
A.l of j37]). Furthermore, a critical point of 5" at fixed 
mass is a (local) maximum iff (see Appendix A.l of [37]): 



S^S=- 



{Spf 
2p 
1 



dO- 



2T 



5pS<Pd9 



mpde] < 0, 



(105) 



for all perturbations Sp that conserve mass: 6M = 0. This 
stability criterion is equivalent to the stability criterion 
P04p but it is simpler because it is expressed in terms of 
the density instead of the distribution function. 

Remark 1: the thermodynamical approach of Kiessling 
|42j in the microcanonical ensemble rigorously justifies the 
maximization problem (J12p . 

Remark 2: comparing the stability criteria (j45p and 
P04p , we see that canonical stability implies microcanon- 
ical stability in the sense that a (local) minimum of F at 
fixed mass is necessarily a (local) maximum of 5' at fixed 
mass and energy. Indeed, if inequality (|45]) is satisfied for 
all perturbations that conserve mass, then inequality p04p 
is satisfied a fortiori for all perturbations that conserve 
mass and energy. However, the reciprocal is wrong in case 
of ensembles inequivalence. Therefore, we just have the 
implication ([T3)) => (HI]). This result can also be obtained 
by comparing the stability criteria (|46l) and (|105p . Indeed, 
since the last term in equation (|105p is negative, it is clear 
that if inequality (j46p is satisfied, then inequality (jlOSp 
is automatically satisfied. In general this is not reciprocal 
and we may have ensembles inequivalence. However, if we 
consider a spatially homogeneous system for which <? is 
uniform, the last term in equation (jlOSp vanishes (since 
the mass is conserved) and the stability criteria (05]) and 
(|105p coincide. Therefore, for spatially homogeneous sys- 
tems, we have ensembles equivalence. 

Remark 3: according to the two interpretations of (J12p 
recalled at the beginning of this section, we note that 
thermodynamical stability implies dynamical stability (for 
isothermal distributions). However, the converse is wrong 
since (|12p provides just a sufficient condition of dynami- 
cal stability. More refined stability criteria are discussed 
in [28129] . 

Remark 4'- the equivalence between (|T^ and (|100p can 
be extended to a larger class of functionals of the form 



S[f] = — ^ Jif^ ^ I) dOdv associated to polytropic dis- 
tributions. Such functionals (Casimirs) arise in the Vlasov 
dynamical stability problem. We refer to |24l43j for a de- 
tailed discussion of this equivalence. 



4.2 The function S{B) 

4.2.1 Global maximization 



The equivalence (|103p is valid for an arbitrary potential of 
interaction u(r, r'). Now, for the HMF model, the problem 
can be simplified further. Indeed, the potential energy is 
given by equation (jH]) so that the energy (|99p and the 
entropy (|102p can be rewritten 



E : 



1,. ttB^ 

-MT — . 

2 k ' 



1 / r2 

S[p] = - I p\npde + -M\n(E+^ 



(106) 



(107) 



We shall first determine the global maximum of entropy 
at fixed mass. To that purpose, we shall reduce the maxi- 
mization problem (jlOOp to an equivalent but simpler max- 
imization problem. 

To solve the maximization problem (|100p . we proceed 
in two steps: we first maximize S[p\ at fixed M and B^ 
and By. Writing the variational problem as 



/ plnpdd \ — aSAI — fi^SB^ — PySBy = 0, 



(108) 



and proceeding as in Sec. 13.2.11 we obtain 

Pi(6i) = Ae-^'=°"(''-'^\ (109) 

where A and A are determined by the constraints Af and 

B through 

M 
A^^^^, (110) 



and 



b = 



27r/o(A) ' 
2ttB /i(A) 



(111) 



kM /o(A)' 

Equation (|109p is the (unique) global maximum of S[p] 
with the previous constraints since S^S = —^ J ^ ^' dO < 
(the constraints are linear in p so that their second order 
variations vanish) . Then, we can express the entropy S as 
a function of B by writing S{B) = S[pi]. After straight- 
forward calculations, we obtain 



SiB)=MlnIo{X)-^-^X + ^ln(E+^ 
k Z \ k 



(112) 



P.H. Chavanis: Dynamical and thermodynamical stability of isothermal distributions in the HMF model 



13 



where X{B) is given by equation piip . Finally, the maxi- 
mization problem (jlOOp is equivalent to the maximization 
problem 



max {5(B)}. 



(113) 



in the sense that the solution of (|100[) is given by equa- 
tions p^ - ([TTT|) where B is the solution of ([TT^ . Note 
that the energy and mass constraints are taken into ac- 
count implicitly in the variational problem (|113p . There- 
fore, (|100[) and ()113p are equivalent for global maximiza- 
tion. However, (|113p is much simpler because, for given E 
and M, we just have to determine the maximum of a func- 
tion S{B) instead of the maximum of a functional S[p] at 
fixed mass and energy. 

Let us therefore study the function S{B) defined by 
equations (J112p and (jllip . Its first derivative is 



S'{B) = M 



/i(A) 27rB 
/o(A) " fcM 



dX 
dB 



2n 



-A 



ttB M 

~E+^ 



(114) 



Using the identity Iq{x) = Ii{x) and the relation (|llip . 
we see that the term in parenthesis vanishes. Then, we get 



^■(«)^?(? 



A 



(115) 



where T is determined by equation (|106p . The critical 
points of S{B), satisfying S'{B) — 0, correspond to 



A 



(3B. 



(116) 



Recalling equations (|106p and (|llip . this leads to the self- 
consistency relations 



2TrB _ hilSB) 
Jm ~ loi^B)'' 



2 k ■ 



(117) 



(118) 



which determine the magnetization as a function of the 
energy. This returns the relationships of Sec. [2l 

Now, a critical point of S{B) is a maximum if S"{B) < 
and a minimum if S"{B) > 0. Differentiating equation 
(jllSp with respect to B, and recalling that the tempera- 
ture T is a function of B given by equation (jl06p , we find 
that 



S"{B) = 



27r / 1 



T 



dX AttB'^ \ 
dB" kMT^ J 



(119) 



Using equation ([59|). we note that 

F"{B) 8t:^B^ 



S"{B) 



Therefore, a critical point of S{B) is a maximum if 

dX 1 inB^ 
dB^ T^ kMT^ ' 



(120) 



(121) 




0.4 0.6 

27tB/kM 

Fig. 8. Entropy S{B) as a function of magnetization for a 
given value of energy. For E > Ec, this curve has a (unique) 
global maximum at 5 = 0. For E < Ec, this curve has a local 
minimum at _B = and a global maximum at B(E) > 0. 



and a minimum if the inequality is reversed. Since the last 
term in equation (|120p is negative, we recover the fact 
that canonical stability implies microcanonical stability. 
Indeed, if the critical point is a minimum of free energy 
{F"{B) > 0), then it is a fortiori a maximum of entropy 
{S"{B) < 0). Furthermore, we know that the series of 
equilibria E{T) is monotonic (sec Sec. [2|). Therefore, the 
homogeneous solution is a maximum of S{B) for E > E^ 
since it is a minimum of F{B) for T > Tc. On the other 
hand, the inhomogeneous solution is a maximum of S{B) 
for E < Ec since it is a minimum of F{B) for T < T^. 
Finally, for the homogeneous solution B = 0, wc have 
S"{B) = -F"{B)/T. Therefore, the homogeneous solu- 
tion is a local minimum of S{B) for E < Ec since it is a 
local maximum of F{B) for T < Tc. Using the equivalence 
between (|100p and (|113p for global maximization, we con- 
clude that the maximum entropy state at fixed mass and 
energy is the homogeneous state for E > Ec and the inho- 
mogeneous state for E < Ec- On the other hand, since the 
ensembles are equivalent for homogeneous solutions (see 
Remark 2 in Sec. 14. ip and since we have established that 
the homogeneous solution is a saddle point of free energy 
for T < Tc (see Sec. I3.2.2p . wc conclude that it is a sad- 
dle point of entropy for E < Ec- This result is shown by 
another method in the next section. 

To complete our analysis, it can be useful to plot the 
function S{B) for prescribed mass and energy. Using equa- 
tions (|112p and (|llip . the normalized entropy s = S/M 
can be expressed in terms of A according to 



s(A)=ln/o(A)-A 






\- 



2h{Xf 

loixr 



(122) 



Eliminating A between expressions (|122p and (|llip . we 
obtain the entropy s{b) as a function of the magnetization 
b for a fixed value of the energy E (more precisely, for 
given e. these equations determine s{b) in a parametric 
form). For E > Ec and E < Ec, this function displays the 
two behaviors described above, as illustrated in Figure [5] 
For E —>■ Ec so that X,b — > 0, we find that the entropy 
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takes the approximate form 



sib) 



1 



In e 



1 



l\b' 



4e2 



-b\ (6^1). 



(123) 



For e < 1 , we explicitly check that the maximum satisfying 
s'{b) = and s"(6) < is given by the second relation in 
equation ([57)) . 

Remark: In Appendix |E1 we plot the second variations 
of entropy S"{B{E)) (related to the variance of the mag- 
netization) as a function of energy and recover the previ- 
ous conditions of stability. 



4.2.2 Local maximization 



We shall now show that the maximization problems (jlOOl) 
and (J113p are also equivalent for local maximization. To 
that purpose, we shall relate the second order variations 
of S[p] and S{B) by using a suitable decomposition. 

A critical point of (|100p is determined by the varia- 
tional principle 



SS - a6M = 0, 



(124) 



where a is a Lagrange multiplier accounting for the con- 
servation of mass. This leads to the distribution (see Ap- 
pendix A. 1. of 127]): 



M 



P 



27r/o(/3B) 



-/3Bcos9 



(125) 



where the temperature is determined by the energy ac- 
cording to equation (|106p . The magnetization B is ob- 
tained by substituting equation (|125p in equations ([S])-® 
leading to the self-consistency relation (P^ . Using equa- 
tions (|105p and ([3]), a critical point of S[p] at fixed mass 
is a (local) maximum iff 



S^S 



M 

2p 



de 



^m.f 



[SByf) 



in^B' 



i5B,f<0, (126) 



for all perturbations 6p that conserve mass: J dpd9 — 0. 
We note that the second order variations of entropy (|126p 
are related to the second order variations of free energy 
13 by 



S^S 



T 



-6'F 



Att^B^ 
k^MT^ 



[SB^f 



(127) 



Writing the perturbation 5p in the form (j68p with equation 
(|5S1) and using expressions (|M)) and ([M|) . we obtain for 
B^Q: 



S^S 



I fF"{B) %f^\ 



de. 



(128) 



and for B ~ 0: 



.2c_ 1^"(0) 



6-'S = --^-^[i5B,r + iSByr 



-If^-^f^de. (129) 

27 p 

Finally, using identity (|120p . we arrive at 

6'S = ls"{B){SB,)' -\l ^^^ de, (130) 
2 2 J p 

for i? 7^ and 

s's = ls"{B) [{SB^r + iSByf] I J ^ de, 

(131) 

for B = 0. These relations imply that p is a local maximum 
of S[p] at fixed mass iff _B is a local maximum of S{B). 
Indeed, if S"{B) < 0, then S^S < since the last term is 
negative. On the other hand, if S"{B) > 0, it suffices to 
consider a perturbation of the form (|68l) with Sp± = and 
Sp\i given by equation ((86)) to conclude that 6^S > for 
this perturbation. This implies that p is not a local max- 
imum of S[p] since there exists a particular perturbation 
that increases the entropy. This is the case for homoge- 
neous solutions when E < Ec- Therefore, ()100p and ()113p 
are equivalent for local maximization. Combining all our 
results, we conclude that the variational problems ()T2)) . 
()100p and ()113p are equivalent for local and global maxi- 
mization 



(IT2])<=» (fTonl) O (fTT3l). 



(132) 



5 Conclusion 



In this paper, we have presented a new method to settle 
the stability of homogeneous and inhomogeneous isother- 
mal distributions in the HMF model. This method starts 
from general variational principles and transforms them 
into equivalent but simpler variational principles until a 
point at which the problem can be easily solved. For 
isothermal distributions, this method returns, as expected, 
the same results as those obtained in the past by different 
procedures |4I29I30J , but we would like to emphasize why 
our approach is interesting and complementary to other 
methods. 

First of all, it is based on general optimization prob- 
lems: the maximization of S[f] at fixed mass M and en- 
ergy E or the minimization of F[f] at fixed mass M. 
These optimization problems provide either conditions of 
thermodynamical stability (in microcanonical and canoni- 
cal ensembles respectively) or sufficient conditions of dy- 
namical stability (more or less refined) with respect to 
the Vlasov equation. Therefore, our approach allows us to 
treat thermodynamical and dynamical stability problems 
with the same formalism. This is not possible if we follow 
an approach starting directly from the density of states 
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g{E) or the partition function Z{/3) which only applies 
to the thermodynamical problem |11I16I30] . As a conse- 
quence, our procedure still works if the Boltzmann func- 
tional is replaced by more general functionals of the form 
S[f] = — J C{f) drdv, where C(/) is convex. In that case, 
the maximization of S at fixed M and E or the minimiza- 
tion of F at fixed M provide conditions of formal dy- 
namical stability for arbitrary steady states of the Vlasov 
equation of the form / = /(e) with /'(e) < [22]. This 
dynamical stability problem has been studied recently for 
polytropic distributions in |24] . Our procedure could also 
be employed for the Lynden-Bell entropy [TSl, although 
the calculations would be more complicated. We have cho- 
sen to treat here the isothermal case in detail since the cal- 
culations are fully analytical. It offers therefore a simple 
illustration of the method. 

Secondly, the optimization problems (fT2|) and (|13|) on 
which our approach is based are both intuitive and rigor- 
ous. The thermodynamical approach is intuitive because 
the Boltzmann entropy S[f] can be obtained from a simple 
combinatorial analysis. In that case, the entropy is pro- 
portional to the logarithm of the disorder where the dis- 
order measures the number of microstates associated with 
a given macrostate. Therefore, maximizing entropy S[f] 
at fixed mass and energy amounts to selecting the most 
probable macroscopic state consistent with the dynami- 
cal constraintfo. In the context of the HMF model, this 
thermodynamical approach has been initiated in |14ll7j 
but was not performed to completion since the stability of 
the inhomogencous phase was not proven (at least analyt- 
ically) by this method. This has been done in the present 
paper. The dynamical approach is also intuitive because 
we qualitatively understand that the stability of a dy- 
namical system is linked to the fact that the system is 
in the minimum of a certain "potential" . For infinite di- 
mensional systems, the nonlinear dynamical stability of a 
steady state of the Vlasov equation relies on the energy- 
Casimir method and its generalizations [52]. In the con- 
text of the HMF model, this approach has been followed 
in |17I25I29| . The optimization problems ^ and ([13]) 
are also rigorous because they have been given a precise 
justification by mathematicians. In statistical mechanics, 
the canonical criterion (|13p has been justified rigorously 
in [TT] and the microcanonical criterion P^ in [32] • En- 
sembles inequivalence has been formalized in [6]. On the 
other hand, the formal and nonlinear dynamical stabil- 
ity of a steady state of the Vlasov equation have been 
discussed extensively in the mathematical literature. We 
refer to [33] for a survey on the standard energy-Casimir 



method and to |45j for refined stability criteria (in the 
context of 2D turbulence) . 

Thirdly, our method determines not only the strict 
caloric curve /?(£') (corresponding to global maxima of S 
at fixed E and M, or global minima of F at fixed M) but 
it also provides the whole series of equilibria containing all 
the critical points of S at fixed E and M, or the critical 
points of F at fixed M. In particular, it allows us to deter- 
nnne metastable states that correspond to local maxima 
of S at fixed E and M, or local minima of F at fixed M. 
There are no such metastable states in the HMF model 
for isothermal distributions, but they can exist in other 
situations |22I24| . Our approach can be used to determine 
the sign of the second order variations of the thermody- 
namical potential in order to settle whether the critical 
point is a global maximum, a local maximum or a saddle 
point. In particular, we have been able to relate the sec- 
ond variations of the different functionals in order to show 
that the equivalence between the optimization problems 
is not only global but also local. 

For these reasons, the approach developed in the 
present paper is an interesting complement to other meth- 
ods [4I16I29I3TJI and it could find application and useful- 
ness in more general situations. 



A Asymptotic expansions 

In this section, we give the asymptotic expansions of b{x)^ 
7/(x), e(a;), s{x) and f{x) for a; — > (corresponding to 
the bifurcation point) and x — > +oo (corresponding to the 
ground state). 

For x — >■ 0, we have 



3 5 



?/(x) 



8 192 



+ o(x6), 



^(-)-l-? + ^+«(-^)' 






(133) 



(134) 



(135) 



(136) 



^ This result can also be derived from kinetic theory. For 
isolated systems, the evolution of the distribution function is 
governed by a kinetic equation that monotonically increases 
entropy S[f] while conserving mass and energy until the maxi- 
mum entropy state is reached. Similarly, for systems in contact 
with a heat bath, the evolution of the distribution function is 
governed by a kinetic equation that monotonically decreases 
free energy F[f] while conserving mass until the minimum free 
energy state is reached [44] . These H-theorems are another way 
to justify the optimization problems (|12|) and ()13p . 



ox 



/(^) = 1-^+0(^6). 



For X — ^ -l-oo, we have 



^(^)- 1-^-8^ + ''^^")' 



(137) 



(138) 



.(-)-fn + Ti + T^ + «(--^)- (139) 
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e(x) = -2+ ^ + o(a; ^), 



X x^ 



(240) C Connection between statistical mechanics 
and variational principles 



s{x)^-\nx + -il-\mT) + j^ 



/(•^) = -2 



2 + 21n7r + 41na; 



In TT + 2 In a; 



+ oix'^). 



/ -3^ ("1 41 "i In this section, we discuss the connection between the sta- 
tistical mechanics of systems with long-range interactions 
and variational principles. 

In the microcanonical ensemble, the accessible config- 
urations (those having the proper value of energy) are 
equiprobable. Therefore, the density probability of the 
(142) configuration (611, wi,...,6'Ar,t;Ar) is P/v(^i, "i, ...,6'Ar, wat) = 
-j^d{E — H) where g{E) is the density of states 



B Derivation of inequality (I104J) 



giE) = / S{E- H)d9idvi...d0Ndv 



N- 



(149) 



We consider a perturbation Sf around a distribution /. 
The exact variations of mass (O and energy ([8|) for any 
perturbation are 



AE 



AM 



6f— dOdv 
■' 2 



5f d0dv, 



(143) 



<PSpde + - I SpS<Pd9dv. 



(144) 



Considering small perturbations Sf, the variations of en- 
tropy (|10p up to second order is 



AS^- /"(hi/ + l)dfd9dv - I 



isfl 

2/ 



d9dv. (145) 



Let us now assume that / is a critical points of entropy at 
fixed mass and energy. It is determined by the variational 
principle ^^, leading to 



\nf + l = -/3 



<P 



(146) 



Substituting this relation in equation (jl45p . we obtain 
„2 X 1 r^5f)2 



AS 



/3 



<? 



5} d9dv 



2/ 



■ d9dv. 



(147) 



Now, using the conservation of mass and energy AAI 
AE = 0, we get 



AS=- 



(Sf) 



1 



d9dv--l3 I 5pS^d9dv. (148) 
Zj z 



Requiring that the critical point be a maximum of entropy 
at fixed mass and energy leads to inequality (|104p . We re- 
fer to [29I46J for generalizations of this result to a larger 
class of functionals in the Vlasov dynamical stability con- 
text. 



In other words, g{E)dE gives the number of microstates 
with energy between E and E+dE. The entropy is defined 
by S{E) = In g{E). Let us introduce the (coarse-grained) 
one-body distribution function f{9, v). A microstate is de- 
termined by the specification of the exact positions and 
velocities {9i,Vi} of the N particles. A macrostate is de- 
termined by the specification of the density {f{9,v)} of 
particles in each cell [9, 9 + d9] x [v, v + dv] irrespectively 
of their precise position in the cell. Let us call f2[f] the 
unconditional number of microstates corresponding to the 
macrostate /. For N — >• -l-oo, equation (|149p can be for- 
mally rewritten 

g{E) ~ / m]5{E[f] - E)S{M[f] - M) Vf, (150) 



where E[f] is the mean field energy ^ and M[f] is the 
mass ([7]). The entropy of the macrostate / is defined by 
the Boltzmann formula S[f] — lni7[/]. The Boltzmann 
entropy can be obtained by a standard combinatorial anal- 
ysis leading, for N -^ +00, to equation ((TU]). Therefore, 
using i7[/] = e'^I^l, equation (|150l) can be rewritten 



liE) 



e^[^l<5(i?[/] - E)diM[f] - M) Vf. (151) 



The unconditional density probability of the distribution 
/ is Po[f] = ^e'^I-'"] (where A is the hypervolume of the 
system in phase space). The microcanonical density prob- 
ability of the distribution / is P[f] = ^e^[/l5(£;[/] - 
E)S{M[f]-M). 

Integrating over the velocities in equation (11491) , a clas- 
sical calculation leads to 



9iE) 



2^ 

r(f) 



[2{E-U)]^^ d9i...d9N, (152) 



where U{9i, ...,9^) is the potential energy (second term 
in the r.h.s. of equation ([T|)). Let us introduce the (coarse- 
grained) one-body density p{9) and denote by f2[p] the 
unconditional number of microstates {9i} corresponding 
to the macrostate p. For N — >■ -l-oo, equation p52p can be 
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formally rcwrittcrPi 

9 

/?(B)=e-"^^+*^'"^°(^). (160) 



The computation of i?(B) is classical [J] and is briefly 

/reproduced in Appendix [DJ with some complements, for 
e"5- '"(^-^[''l)f2[p],5(A/[p] - M) Vp, (153) the sake of self-consistency. For N -^ +cx), we have 



where W[p] is the mean field potential energy (second 

term in the r.h.s. of equation (|5|)). The unconditional 

number of microstates f2[p] can be obtained by a clas- Therefore, the density of states (|159[) can finally be written 

sical combinatorial analysis leading, for N — >■ +oo, to 

n[p] = g-/pinpde^ Therefore, equation p35|) can be g(-B) ~ /"e^^-^^dB, (161) 

rewritten J 



giE)^ eS^P^d{M[p]-M)Vp, 



(154) where S{B) is giyen by equation (J112p . The unconditional 
density probability of the magnetization B is Po(B) = 



inhere S[p] is giyen by equation (|102p . The uncondi- 



-A+Min/o(A)^ The microcanonical density probabil- 



tional density probability of the distribution p is Po[p] = ity of the magnetization B is P(B) = -t^g'^^-^'. 

^g-J pinpde ^^Yiere A is the yolume of the system in In the canonical ensemble, the density prob- 



physical space). The microcanonical density probability ability of the configuration 



h,vi, 



yN 



vn) is 



of the distribution p is P[p] 



JiE) 



zSM5{M[p] - M). 



PN{di,vi,...,6N,VN) = zTm^ '^^ where Z(/3) is the 



For the HMF model, let us introduce the magnetiza- partition function 
tion yector 



z{P) 



B-r ^ —- 



k 



N 



N 



m-Je-^" 



d9idvi...d9Ndv 



•n- 



. ,,.... (162) 

■V^cos^i, By = ---y^j,\nBi. (155) ■' 

*~ '~ The free energy is defined by F{P) = — 4lnZ(/3). Intro- 

Z!^!Z?*„''^!*/„''.^ 1''''"^^ '''''' ^"^ expressed in terms of the ^^^j^g ^^^ unconditional number of microstates Q\f] cor- 

responding to the macrostate /, we obtain for N — > -|-cx): 



magnetization as 



U 



ttB^ 



kN 
47r ' 



(156) 



Z{(3) 



Therefore, the density of states (|152p can be rewritten 



.-0Elf] 



n[f]S{M[f]-M)Vf 



g{E) 



27r^ 



r{§ 



2 E 



ttB^ 



kN\ 
4^ J 



eS[f]~mf]s{M[f]-M)Vf 



e-l'^P^fU{M[f]-M)Vf, (163) 



xS 



(2ttB^ 






where F[f] ^ E[f] — TS[f] is the Boltzmann free energy 
(|lip . The canonical density probability of the distribution 

/isP[/] = ^e-^^[/l<5(M[/]-M). 

Integrating oyer the yelocities in equation (|162p . we 
(157) get 



or, equiyalently, 

2^^/2 



9{E)^ 



r(f) 
/27rS^ 2ttB. 
\ k '~r 



2 E 



ttB' 



kN 
An 



m-i'i 



N/2 



e-P'^ d9i...d0 



N- 



(164) 



<^)<^ 



Introducing the unconditional number of microstates f2[p] 
corresponding to the macrostate p, we obtain for A^ — > 

+oo: 



(158) 

where ^(B) denotes the unconditional number of mi- 
crostates {9i} corresponding to the macrostate B. For 
N — >■ -l-oo, we haye 



9iE) 



gf ln(S+^ 



'^{B)dB. 



(159) 



Z{I3) ~ / e-'^^^P^n{p]SiM[p] - M)Vp 
3- 1 "''"""''^^^''UiMlp] ~ AI)Vp 

e-0P^PU{M[p]-M)Vp, (165) 



Like in the main part of the paper, for the sake of concise- 
ness, we do not explicitly write the constant terms that are 
independent on p in the expression of the entropy (the term in 
the exponential). They can be restored easily. 



where F[p] is given by equation P5|) . The canonical 
density probability of the distribution p is P[p\ = 



zm^ 



-I^Flp]S{M[f] - A/). 
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For the HMF model for which the potential energy can 
be expressed in terms of the magnetization, the partition 
function (|164p can be rewritten 



-^(l)"'7 



-I3U 




smf 



.,(!!^),(M.|,,,^^,,,„, 



(166) 



or, equivalent ly. 



^(«H'f 



•xQ 



2ttB^ 2ttB. 



^)<^) 



'27rB, 



k ■ k ,-~y , /"["t^^- ""> 



For N — > +00, we get 

Z(^)~ f e^'^f2{B)dB 



(168) 



where F{B) is given by equation ([55)) . The canonical 
density probability of the magnetization B is P(B) = 

_J_f,-l3F{B) 

Let us now denote by </> a generic global variable such 
as f{6,v), p{9) or B. We also recall that for systems 
with long-range interactions, for which the mean field ap- 
proximation is exact in the proper thermodynamic limit 
N — !• -l-oo, we have the extensive scalings S[cl)] = Ns[(j)], 
E[(j)] = Ne[(j)], F[(j)] = Nf[(j)]. Accordingly, the preceding 
results can be formally written 



and 



9iE) 



Z{P) 



-,Ns\, 



^5{cmce[4>]~cmce)V(J>, (169) 



^-PNm si^ccE^ - cce) ■Dcj), (170) 



where the (5-functions take into account the con- 
straints as described above. The microcanonical num- 
ber of microstates corresponding to the macrostate (f> 

is i7[0] ~ e^^I*^! 5(cAfC£;[0] — cmce) and the micro- 
canonical probability of the macrostate (j) is P[4'\ = 
-^e^^l"^! J(ca/cb[0] — cmce)- Similarly, the canonical 
number of microstates corresponding to the macrostate 
(j) is n[(l)] = e-'^^^[*U(cc£M - cce) and the 
canonical probability of the macrostate (j> is ^[0] = 
.e-/3^/[0] S{cce[<I>]-cce)- For N -^ +00, we can make 



1 



the saddle point approximation. In the microcanonical en- 
semble, we obtain 



9{E) 



.SiE) _ ^Nsl,p,] 



I.e. 



lini ^SiE)^s[^4, 



(171) 



(172) 



where 0* is the solution of the maximization problem 

mayi{s[(l)]\E,M}. (173) 



This leads to the variational problems (J12p . (|100[) and 
(|113p . In the canonical ensemble, we obtain 



Z(/3)=e-^^(^)~e-'3^^[^*l, 



hm 1f(/3)=/[0,] 



(174) 



(175) 



where 0* is the solution of the minimization problem 

min{/[0] I M}. (176) 



This leads to the variational problems (|T3)) . (|42|) and ([54]) . 
The preceding discussion shows the connection be- 
tween the statistical mechanics of systems with long-range 
interactions (based on the calculation of the density of 
states and of the partition function) and variational prin- 
ciples (based on the maximization of entropy or minimiza- 
tion of free energy). It also shows how the variational 
problems ^, (jTUD)) . ([TT^ and ^, (gl]), ^ are re- 
lated to each other. These results can be made rigorous 
by using the theory of large deviations. We refer to Ellis 
[TT] for a mathematical presentation of this theory and 
to Barre et al. [301 and Touchette [48| for its application 
to physical problems. In the present paper, we have con- 
sidered a different approach. Wc started from the funda- 
mental variational problems ([T^ and (|13p that can be 
motivated by a simple combinatorial analysis. This is the 
historical approach of the problem finding its roots in 
Boltzmann's work. This is also the traditional approach 
used by physicists to determine the statistical equilibrium 
state of self-gravitating systems (see, e.g., [8110149150] ). 
two-dimensional point vortices (see, e.g., jllSlj ) and the 
HMF model (see, e.g., |14I17) ). Therefore, describing the 
statistical mechanics of systems with long-range interac- 
tions from the fundamental variational problems ((T^ and 
P3p and reducing them to simpler but equivalent forms 
([TUP)) . pT^ . P^ and dm) as we have done here is an in- 
teresting presentation that complements the one followed 
in j3]. A bonus of this approach is that it remains valid 
when the variational problems p^ and ([T^ have a dy- 
namical interpretation in relation to the (formal) nonlin- 
ear stability of the system with respect to the Vlasov equa- 
tion [55]. In that case, S\f] is a Casimir functional of the 
form S[f] = — J C{f) dOdv (sometimes called a pseudo 
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entropy) that is more general than the Boltzmann fune- 
tional arising in the thermodynamieal approaeh. In par- 
tieular, in the dynamieal stability problem, the variational 
problems P^ and (fT5|) cannot be obtained from a theory 
of large deviations since their physical interpretations are 
completely different. 

Remark 1: the density of states (|149|) and the partition 
function (|153p can also be calculated from field theoret- 
ical methods (see e.g. Horwitz & Katz [55] and de Vega 
& Sanchez |53] for self-gravitating systems and Antoni & 
Ruffo [TS] for the HMF model). These approaches are valu- 
able but they are also considerably more abstract than the 
one based on variational principles. 

Remark 2: the mean field Boltzmann distribution (|16p 
can also be obtained from the first equation of the Yvon- 
Bogoliubov-Green (YBG) hierarchy [^, but the condition 
of stability (related to the correlation functions appearing 
in the next equations of the YBG hierarchy) is more diffi- 
cult to obtain than with the approach based on variational 
principles. 

Remark 3: for N — > +oo, the density of states (resp. 
partition function) is dominated by the global maximum 
of entropy at fixed mass and energy (resp. minimum of 
free energy at fixed mass) according to (|17ip (resp. (|174[) '). 
Nevertheless, ZocaZ entropy maxima (resp. free energy min- 
ima), i.e. metastable states, are also fully relevant because 
they have very long lifetimes scaling like e^ |55l56j . This 
is particularly true in the case of classical self-gravitating 
systems for which there is no global entropy maximum 
(resp. free energy minimum) [8110] . 



D Distribution of the magnetization 

In this Appendix, we determine the distribution of the 
magnetization by a direct calculation and show its con- 
nection with the entropy S{B) and free energy F{B). The 
density probability of the magnetization vector defined by 
equation (|155p is 



'27^8^ 27rR 



X(5 



'2ttB,, 



= S 



2ttB,: 



^ COS Oi I 



E 



smf 



P, 



N[ 



9N)ddi...deN. (177) 



In the microcanonical ensemble, the A'^-body distribution 
function is Pn{9i,vi, ...,9n,vn) = -tWi5(£' — H). Inte- 
grating over the velocities, a classical calculation gives 



P, 



N[ 



^N 



1 2^^/^ 

9iE)rm 



[2{E-U)]^, (178) 



where U is the potential energy (second term in the 
r.h.s. of equation ([IJ). In the canonical ensemble, the 
A^-body distribution function is P/v(6'i, wi, ..., 6'jv, wat) = 
_i g-^ff ^ Integrating over the velocities, we obtain 



Recalling that the potential energy can be expressed in 
terms of the magnetization according to equation (jl56p . 
the density probability of the magnetization in micro- 
canonical and canonical ensembles is given by 



p,'27:B^ 2Tr By 



1 2^^/2 



2 E^ 



ttB' 



kN 
An 



k J ' 9{E)r{^) 
(2ttB^ 2Tr By' 



Q 



V k 



(180) 



p ^ 2nB, 2TTBy '^ ^ _J_ f2TT_'^ 



N/2 



\ k ' k J Z{P) \j3 J 
p{^^^)^(2nB,, 2nBy 



xe 



^Q 



\ k ' k 



(181) 



where 




xJ I lll^+^sin^i j d6i...de_ 



N, 



(182) 



is the unconditional number of microstates with magneti- 
zation B. The calculation of this integral is classical [3]. 
Using the Fourier representation of the i5-function 



'"^°° da 
Zn 



(183) 



we obtain 



n = (27r)^-2 / dq^dg^e^''(9..9.)^ (134) 



where 



2tt 
h{Qx,qy) = i — {qxBj, + qyBy) + \n Jo{q), (185) 



and q = wg2 + g2 Recalling that k ~ 1/A^, the function 

h docs not depend on A^. For A^ — > -f-oo. we can make the 
saddle point approximation 



{] ^ eNh{q:,q;) 



(186) 



where {q*,q*) corresponds to the maximum of h[qx,qy). 
The vanishing of dh/dqx and dh/dqy leads to 






(187) 



P, 



N[ 



yw) = 



Z{P) 



N/2 



-pu 



(179) 



27r Ji{q) qy _ 

^^By-^^-^0, 



Mq) q 



(188) 
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where we have used Jq{x) = — Ji(x). Setting q^ 
and Qy = —i\y, we find that q ~ iX where A = 



-iA.r 



Arr -\- Ay. 



Substituting these expressions in equations (|187p and 

/o(A), we 



(189) 



(jlSSp and using Ji(iA) = i/i(A) and Jo(JA) 
find that A is determined by 



27rB _ /i(A) 
fcM ~ /o(A)' 



Then, we get 



E Variance of the magnetization 

In this Appendix, we show the connection between the 
variance of the magnetization in canonical and micro- 
canonical ensembles and the second order derivatives of 
entropy and free energy. 

The unconditional density probability of the magneti- 
zation (|192p can be written 



^o(B) 



A 



(197) 



<lx 









(190) 



Substituting these values in equations (|185p and (|186p . we 
finally obtain 



/?(B) 



o^^hw^+ln-ToCA)] 



(191) 



Therefore, for N — >■ -|-c», the unconditional density prob- 
ability of the magnetization, i.e. the one corresponding to 
a Poissonian (uncorrelated) distribution of angles, is 



where we have introduced the function 



h{B) 



27rB 



A + ln/o(A), 



(198) 



where \{B) is defined by equation (|189p . Note that h{B) 
can be interpreted as the entropy of the magnetization for 
a Poissonian distribution of angles. For N — >■ -|-oo, the 
distribution is strongly peaked around its maximum. The 
most probable value of the magnetization B corresponds 
to the maximum of h{B). Using equation (|189p . the first 
derivative of h{B) is 



A(B) - ^e 



N[-'W->'+^^I«W] 



(192) 






(199) 



This result can also be obtained from the theory of large 
deviations by a direct application of the Cramer theorem 
(see equations (6) and (7) in [50]). 

We now take into account the correlations by using 
the A^-body distribution functions ([T751) and ^75)) . Ac- 
cording to equations (jl80p and ()19ip the distribution of 
the magnetization in the microcanonical ensemble is given, 
for N — >• -l-oo, by 



1 ^^ In ( t^j^jiMjL 



(£;+^)-2^A+Mln7o(A) 



It can be written 



Pmce{B) 



9{E) 



^S{B) 



(193) 



(194) 



where S{B) is the entropy defined in equation pi2p . Ac- 
cording to equations (|18ip and (J19ip . the distribution of 
the magnetization in the canonical ensemble is given, for 

A^ -J> -l-oo, by 



Pcb(B) 



1 gi^-^A+MIn/„(A)^ (^95) 



zm 



It can be written 



PCE^B) 



Z{P) 



'I3F{B) 



(196) 



where F{B) is the free energy defined in equation ([5 



so that the most probable value of the magnetization is 
A = B = (we will see that it corresponds indeed to 
a maximum). Let us expand h[B) around its maximum 
/i(0) — 0. The second derivative of h{B) is 



h"{B) = - 



27r d\ 
kMdB' 



Using identity 



we obtain 



h"iO) =. - 



8n' 



<0, 



(200) 



(201) 



justifying that _B = really is the maximum of h{B). 
Therefore, for ^/NB ^ 1, the distribution of the magneti- 
zation is the Gaussian: 



^o(B) 



An 






(202) 



This result can be directly obtained from the central limit 
theorem. The variance of the unconditional distribution 
of magnetization is 



{B')o = 



k^M 

47r2 ■ 



(203) 



Let us now consider the distribution of the magneti- 
zation in the canonical ensemble given by equation (jl96p . 
The most probable value of B corresponds to the min- 
imum of free energy F{B) as studied in section 13.2.11 
The vanishing of F'{B) leads to equation ((58| . The sec- 
ond derivative of F{B) at the extremum point is given by 
equation ([5^ where X'{B) is given by equation (|5^ . 
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In the homogeneous phase -6 = 0, the variance of the 
magnetization is 



(B') 



CE 



(3F"{0)' 
Using equation (j82p . we obtain 

47r 

According to equation ([5^. we have 



(204) 



/3^"(0) 



8n' 



Tc 

1 

T 



Therefore, using equation (|204p . we obtain 



{B') 



kHi 



CE 



47r2 1 - Tc/T ' 



(206) 



(207) 



We first note that, according to equation (j206p . the homo- 
geneous phase is a minimum of free energy (_F"(0) > 0) 
for T > Tc and a maximum of free energy (F"(0) < 0) 
for T < Tc in agreement with the graphical construction 
of section 13.2.11 On the other hand, for T -> -|-cxj, i.e. 
/3 = 0, the distribution of angles becomes uniformly dis- 
tributed and uncorrelated (see equation (|179p ) so that we 
recover the result (|203p valid for a Poissonian distribution. 
Finally, we note that the variance diverges for T — >■ T^ . 
This result was previously obtained in |17l54j (expressed 
in terms of the variance of the force {F'^) = {B'^)/2) from 
the second equation of the YBG hierarchy. 

In the inhomogeneous phase S ^ 0, the variance of 
the magnetization is 



{{ABf 



1 



ICE 



PF"{B) 

Using equations (|82|) and ((58)) . we obtain 

1 _ kM 2nB^ 

X'{B) " "2^ " " kM ■ 

According to equation ([5^. we have 

/3F"{B) 



BM \2-Sr 






(208) 



(209) 



(210) 



Therefore, the variance of the magnetization is given by 
equations (|208p and (j210p where B{T) is given by equation 
((58)) . Figure [9] shows that F"{B{T)) is always positive so 
that the inhomogeneous phase is always a minimum of 
free energy in agreement with the graphical construction 
of section 13.2.11 

Let us finally consider the distribution of the magneti- 
zation in the microcanonical ensemble given by equation 
((194p . The most probable value of B corresponds to the 
maximum of entropy S{B) as studied in section 11.2. II The 
vanishing of S'{B) leads to equation ()117p with equation 



(|118p . The second derivative of S{B) is related to the sec- 
ond derivative of free energy by equation ((120p . 

In the homogeneous phase i? = 0, the variance of the 
magnetization is 



(B'] 



MCE 



S"iO) 



(205) Using equations ((120p and ((206p . we obtain 



S'\Q) = -^^"(0) = - 



_8^ 
k^M 



1 - 



Tc 
T 



(211) 



(212) 



Using equation pisp with 5 = 0, this can be expressed 
in terms of the energy as 



S"iO) 



Stt^ / Ec 
k^M \ E 



Therefore, using equation ((21ip . we obtain 



iB' 



k^M 



>MCE 



47r2 \~EclE' 



(213) 



(214) 



We note that, according to equation (|213p . the homo- 
geneous phase is a maximum of entropy (S"'(0) < 0) 
for E > Ec and a minimum of entropy (S"'(0) > 0) 
for E < Ec in agreement with the discussion of section 
14.2.11 On the other hand, the variances of the magnetiza- 
tion in canonical and microcanonical ensembles coincide: 
{B )cE = {B )mce- 

In the inhomogeneous phase B ^ Q, the variance of 
the magnetization is 



{{ABf 



'MCE 



S"{E)' 
Using equations (|120p and ()210p . we obtain 



(215) 



S"{B) = - 



87r2 
BM 



1 



9 _ X _ SiT^B- 




(216) 



Therefore, the variance of the magnetization is given by 
equations ()215p . ()216p where B{T) is given by equation 
(|117p . It can be expressed in terms of the energy by using 
equation ()118p . Figure fTOJ shows that S"{B) is always neg- 
ative so that the inhomogeneous phase is always a maxi- 
mum of entropy in agreement with the results of section 
14.2.11 The variances of the magnetization in canonical and 
microcanonical ensembles do not coincide in the inhomo- 
geneous phase: {{ABY)ce t^ {{ABY)mce- 

In order to represent these results graphically (see Fig- 
ures and nop , it is convenient to introduce the dimen- 
sionless variables defined in section |5| The variance of the 
magnetization in the homogeneous phase is 



N{b^ 



1 



1 



1-1/e' 



(217) 
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Fig. 9. Variance of the magnetization as a function of the 
inverse temperature. We have represented the variance in the 
canonical (full line) and microcanonical (dashed line) ensem- 
bles. 



On the other hand, close to the ground state (e, 77) — > 

(-2+,+(»): 



N{{Ahf)cE 



N{{Ahf 



>CE 



1 



16772 



32 



(g + 2)^ 
64 



(222) 



(223) 



We emphasize that althought the ensembles are equiv- 
alent regarding the caloric curve P{E), the variance of 
the magnetization in the inhomogeneous phase differs in 
the two ensembles. Therefore, numerical simulations of the 
isolated HMF model (microcanonical ensemble) and of the 
dissipative BMF model (canonical ensemble) should lead 
to different values of {{AB)^). 



10 r 



^ 4- 



CE / 



/MCE 



Fig. 10. Variance of the magnetization as a function of the 
energy. We have represented the variance in the microcanonical 
(full line) and canonical (dashed line) ensembles. 



both in canonical and microcanonical ensembles. The vari- 
ance of the magnetization in the inhomogeneous phase is 



N{{Abf)cE 



1 



2ri 



2rj' 



in the canonical ensemble and 



N{{Abr] 



MCE 



j^^^ - 277 + 8b^v- 



(218) 



:, (219) 



in the microcanonical ensemble. In the microcanonical en- 
semble, it can be expressed in terms of the energy, us- 
ing e = 1/77 — 2b^. Using the asymptotic expansions of 
Appendix |21 we obtain close to the bifurcation point 

(e,77)^(l-,l+): 



N{iAbf)cE 



N{{Abf)MCE 



1 



4(r;-l) 4(1 -e)' 



20(?y-l) 4(1 -e)' 



(220) 



(221) 
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